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Abstract 



The main arguments in the discussion of the proper treatment of the total inclusive 
cross section for bottom-quark associated Higgs-boson production are briefly reviewed. 
A simple and pragmatic formula for the combination of the so-called four- and five- 
' flavour schemes is suggested, including the treatment of the respective theory error 

, estimates. The numerical effects of this matching formula are discussed. 

rn 

^ ; 1 Two approaches 

The inclusive total cross section for bottom-quark associated Higgs-boson production, de- 
note4l] pp/pp {bb)H + X, can be calculated in two different schemes. As the mass 
^ , of the bottom-quark is large compared to the QCD scale, mb ^ AqcD) bottom-quark 

^ I production is a perturbative process and can be calculated order by order. Thus, in a 

four-flavour scheme (4FS), where one does not consider 6-quarks as partons in the proton, 
the lowest-order QCD production processes are gluon-gluon fusion and quark-antiquark an- 
nihilation, gg bbH and qq — > bbH, respectively. However, the inclusive cross section for 
gg {bb)H develops logarithms of the form \n[^-p/m\,), which arise from the splitting of 
gluons into nearly collinear bb pairs. The large factorization scale /ip ~ rn^i/A corresponds 
to the upper limit of the collinear region up to which factorization is valid [U [21 [3] . For 
Higgs-boson masses ttt-h ^ 4mb, the logarithms become large and spoil the convergence 
of the perturbative series. The \n.{pi-p / m^,) terms can be summed to all orders in perturba- 
tion theory by introducing bottom parton densities. This defines the so-called five-flavour 
scheme (5FS). The use of bottom distribution functions is based on the approximation 
that the outgoing 6-quarks are at small transverse momentum. In this scheme, the LO 
process for the inclusive {bb)H cross section is bottom fusion, bb — ?• H. 

If all orders in perturbation theory were taken into account, the four- and five-flavour 
schemes would be identical, but the way of ordering the perturbative expansion is different. 



^The notation {bb)H is meant to indicate that the 66 pair is not required as part of the signature in 
this process, so that its final state momenta must be integrated over the full phase space. 
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At any finite order, the two schemes include different parts of the all-order result, and the 
cross section predictions do thus not match exactly. While this leads to an ambiguity in 
the way the cross section is calculated, it also offers an opportunity to test the importance 
of various higher-order terms and the reliability of the theoretical prediction. The 4FS 
calculation is available at NLO [HE], while the 5FS cross section has been calculated at 
NNLO accuracy [6]. Electroweak corrections to the 5FS process have been found to be 
small [7] and will not be considered in the numerical results presented in this note. 

In the next section we briefly summarize some of the features of the two schemes. 

2 Discussion 

In this section, we collect the main arguments and counter-arguments that have been 
discussed in the physics community for and against either the 4FS or the 5FS approach. Let 
us stress that at this point, we consider only the total inclusive cross section pp — > {bb)H + 
X. Once distributions or even tagged ^-quarks are taken into account, the discussion below 
has to be re-evaluated. 

Gluon to b-quark splitting — kinematical approximations and collinear logarithms: In pp 
or pp collisions, bottom-quarks arise at LO in through the splitting of a gluon. In 
the 4FS, this splitting is described by fixed-order perturbation theory, and includes 
the full dependence on the transverse momentum px of the 6-quark and its mass. In 
the 5FS scheme, the splitting arises by solving the DGLAP evolution equations with 
five massless quark flavours [8]. 

Integrating over phase space, the 5FS at LO neglects all contributions beyond 
as pt — >• 0, while they are consistently taken into account in the 4FS. However, at 
higher orders, the 5FS does include sub- leading flnite pT effects; at NNLO, all the 
information of the LO 4FS calculation is incorporated in the 5FS as well. Note, 
however, that the 4FS result is available through NLO, so that the finite px effects 
are available one order higher than in the 5FS calculation. 

Since the factorization theorem holds only for massless quarks, in the 5FS scheme 
one has to assume mb = (the bottom Yukawa coupling can - and must - be kept 
finite though). As the only other scales of the inclusive cross section are mn and the 
hadronic center-of-mass energy s ^ m^, the 4FS and the 5FS formally differ at order 
m^/m^, even if all orders of the perturbation theory could be taken into account. 
The numerical size of these effects has been investigated in Ref. [9] and was found 
to be negligible. 

In the opposite phase space region of very small transverse momenta of the 6-quarks, 
the 5FS implicitly re-sums logarithmic terms of order a^^"^ /p^ ■ \v?^~^ p^ (n > 0), 
where A; = in the LO, k G {0,1} in the NLO, and k € {0,1,2} in the NNLO 
calculation. In the 4FS, these effects are not re-summed, but taken into account 
only up to finite n. 

Renormalization/ factorization scale dependence: Due to the different complexity of the 
two LO processes (2^3 for the 4FS and 2—7-1 for the 5FS), radiative corrections 
are known through NNLO in the 5FS [6], while they are available only through NLO 
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in the 4FS [HIS]. This results in a stronger dependence in the 4FS on the unphysical 
renormahzation and factorization scales, /xr and fip, than in the 5FS. 

Top-quark induced effects: Currently, the 5FS calculation neglects effects proportional to 
the top Yukawa coupling. Partly, they are taken into account in the gluon fusion 
induced cross section, gg ^ H + X . What is missing are interference effects pro- 
portional to both the top and the bottom Yukawa coupling. In most models with 
an enhanced bottom Yukawa coupling (like SUSY at large tan/3), the top Yukawa 
coupling is simultaneously suppressed though, and therefore these interference terms 
can be neglected. For consistency, these terms have been removed also from the 4FS 
calculations in the numerical results below. 

Parton distribution functions: Technically, the transition from Uf to uj + 1 flavours is 
well understood through NNLO [TU]. The exact value /Uthr where this transition is 
imposed is arbitrary (but fixed in all available PDF sets). At fixed order in perturba- 
tion theory, this introduces an uncertainty in the bottom-quark PDFs which should, 
however, be small for mn ~ //f ^ A'thr ~ f^h- 

In fact, until recently, modern PDF sets always assumed a five- flavour content of 
the proton. This led to an inconsistency in the 4FS calculation. With the return of 
four-flavour PDF sets [111 112^ [T3]. this inconsistency could be removed. 

Numerical phase space integration: Phase space integration in the 4FS calculation is done 
numerically. In the collinear region, this requires high precision and is computation- 
ally very demanding. In the 5FS, the analytical result for the phase space integration 
is published and implemented in a publicly available program bbhSnnlo ^ which 
performs the convolution with the PDFs. The calculation of the cross section for a 
single set of parameters takes of the order of seconds with IbbhOnnlol 

Disagreement between the 4FS and 5FS: If the total inclusive cross section is calculated at 
LO in both the 4FS and 5FS, and the factorization scale is set to /xp = mn, one 
finds that the numerical results disagree by a factor of the order of five (see, e.g., 
Ref. [H]). As already mentioned in Section [H it was argued that in the 5FS, the 
"proper" choice of the factorization scale is /ip ~ '"T'h/4. In fact, the 4FS and 5FS 
results agree very well in this case. The NNLO calculation in the 5FS shows a rather 
mild dependence on /ip, and indeed confirms that /xp = mH/4 is a reasonable choice, 
since the radiative corrections exhibit a good convergence behaviour for that value. 

3 Santander matching 

Theoretically, the 4FS and the 5FS are equivalent descriptions of the total inclusive cross 
section for (bb)H production. However, due to the arguments listed above, differences in 
the numerical predictions are observed at finite order of perturbation theory. By making 
what we believe are reasonable assumptions, in this section we suggest a way to combine 
both approaches. 

The 4FS and 5FS calculations provide the unique description of the cross section in the 
asymptotic limits mn/m]^ — > 1 and mn/jTib oo, respectively. For phenomenologically 
relevant Higgs-boson masses away from these asymptotic regions both schemes are appli- 
cable and include different types of higher-order contributions. The matching we suggest 
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below interpolates between the asymptotic limits of very light and very heavy Higgs- 
bosons. It is pragmatic and far from being theoretically rigorous. However, we believe 
that the matching catches the essence of the arguments given in Section [2] and provides a 
prediction with an uncertainty band which covers the physical result. 

A comparison of the 4FS and 5FS calculations reveals that both are in numerical agree- 
ment for moderate Higgs-boson masses (see Fig. 23 of Ref. [E])- Once larger Higgs-boson 
masses are considered, the effect of the collinear logarithms ln(mH/?7ib) becomes more and 
more important and the two approaches begin to differ. We suggest to combine the two 
approaches in such a way that they are given variable weight, depending on the value of 
the Higgs-boson mass. 

The difference between the two approaches is formally logarithmic. Therefore, the de- 
pendence of their relative importance on the Higgs-boson mass should be controlled by a 
logarithmic term. We determine the coefficients such that 

(a) the 5FS gets 100% weight in the limit mn/mh oo ; 

(b) the 4FS gets 100% weight in the limit where the logarithms are "small". There is 
obviously quite some arbitrariness in this statement. We assume here that "small" 
means ln(?nH/w,b) = 2. The consequence of this particular choice is that the 4FS and 
the 5FS both get the same weight for Higgs-boson masses around 100 GeV, consistent 
with the observed agreement between the 4FS and the 5FS in this regionJl 

This leads to the following formulcH 

matched _ ^ W (T 

l + W 



with the weight w defined as 



1 mn 

w = m 

mb 



and o"^^^ and a^^^ denote the total inclusive cross section in the 4FS and the 5FS, respec- 
tively. For ?7ib = 4.75 GeV and specific values of mn, this leads to 

^matched I n /in^4FS , n „5FS 

^ Lh=iooGcV = 0-^9c7 +0.51C7 , 

^matched | g gg ^4FS ^ o.64 CJ^^^ , 

I mn =200 GeV ' 

^matched | g ^4FS ^ o.69 CJ^F^ , (3) 

I mn =300 Gov ' 

^matched I n nn. „4FS , n Vl ^5FS 

I mn =400 GeV ' 
^matched I n ov ^4FS , n 11 „5FS 

Knnr- -Sf = 0.2/CJ +0.(60 

I mn =500 GeV 

A graphical representation of the weight factor w is shown in Fig.[T](a). 

Concerning the uncertainties, we suggest to add them linearly, using the weights w defined 
in Eq. ([2|) . This ensures that the combined error is always larger than the minimum of the 



^Note that one should use the pole mass for mb here rather than the running mass, since it is really 
the dynamical mass that rules the re-summed logarithms. 

^This formula originated from discussions among the authors at the Higgs Days at Santander 2009 and 
is therefore dubbed "Santander matching". 
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two individual errors. Neglecting correlations and assuming equality of the uncertainties 
in the 4FS and 5FS calculations would imply that the matched uncertainty is reduced by 
a factor of w/{l + w) with respect to the common individual uncertainties. This seems 
unreasonable. In our approach the matched uncertainty would be equal to the individual 
ones in this case. On the other hand, taking the envelope of the 4FS and 5FS error bands 
seems overly conservative to us. 

The theoretical uncertainties in the 4FS and the 5FS calculations are obtained through 
Hp-, PDF-, and variation as described in Ref. [TH]. They can be quite asymmetric, 
which is why the combination should be done separately for the upper and the lower 
uncertainty limits: 

Aa± = ^ (4) 

I + w 

where Ac^f^ and Ac^F^ are the upper /lower uncertainty limits of the 4FS and the 5FS, 
respectively. 



4 Numerical results 

We shall discuss the numerical implications of the matching prescriptions defined in the 
previous section and provide matched predictions for the inclusive cross section pp — >■ 
{bb)H + X at the LHC operating at a center-of-mass energy of 7 TeV. The individual 
numerical results for the 4FS and 5FS calculations have been obtained in the context 



of the LHC Higgs Cross Section Working Group with input parameters as described in 



Ref. [E]. Note that the cross section predictions presented below correspond to Standard 
Model bottom Yukawa couplings and a bottom-quark mass of rrib = 4.75 GeV. SUSY 
effects can be taken into account by simply rescaling the bottom Yukawa coupling to the 
proper value [71 [16] . 

Fig.[T](b) shows the central values for the 4FS and the 5FS cross section, as well as the 
matched result, as a function of the Higgs-boson mass. The ratio of the central 4FS and 
the 5FS predictions to the matched result is displayed in Fig.[2](b) (central dashed and 
dotted line): for mn = 100 GeV, the 4FS and the 5FS contribute with approximately the 
same weight to the matched cross section, with deviations between the individual 4FS and 
the 5FS predictions and the matched one of less than 5%. With increasing Higgs-boson 
mass, the 4FS result deviates more and more from the matched cross section due to its 
decreasing weight. At mu = 500 GeV, it agrees to less than 20% with the matched result, 
while the 5FS is still within 8% of the latter. 

The corresponding theory error estimates are shown in Fig.[2j The absolute numbers are 
displayed in panel (a), while in panel (b) they are shown relative to the central value of the 
matched result. Up to mn ~ 300 GeV, the combined uncertainty band covers the central 
values of both the 4FS and the 5FS. For larger Higgs-boson masses, the 4FS central value 
is slightly outside this band. 
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Figure 1: (a) Weight factor w, Eq. ([2]), as a function of the Higgs-boson mass 
mu- The bottom-quark pole mass has been set to mb = 4.75 GeV. (b) Central 
values for the total inclusive cross section in the 4FS (red, dashed), the 5FS (green, 
dotted), and for the matched cross section (blue, solid). Here and in the following 
we use the MSTW2008PDF set [ll] (NLO for 4FS, NNLO for 5FS). 




Figure 2: (a) Theory uncertainty bands for the total inclusive cross section in 
the 4FS (red, dashed), the 5FS (green, dotted), and for the matched cross section 
(blue, solid), (b) Uncertainty bands and central values, relative to the central 
value of the matched result (same line coding as panel (a)). 
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5 Conclusions 



We have reviewed the main arguments for and against following either the 4FS or the 5FS 
when predicting the value of the total inclusive cross section for bottom-quark associated 
Higgs-boson production. Based on this discussion, we have suggested a way to combine 
the two approaches such that their central values and uncertainties enter with appropriate 
weights. 
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